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We study the equation of state of quark matter at finite 
temperature, using a confinement model in which chiral sym- 
metry remains broken in the deconfined phase. Implications 
for type II supernova explosion and for the structure and evo- 
lution of the proto-neutron star are discussed. 

PACS numbers: 26.50,+x, 24.85.+P, 12.39.-x, 26.60.+C 



The Equation of State (EOS) of matter at high den- 
sity and temperature is a very interesting but contro- 
versial topic. Among the possible applications of such an 
EOS are the structure of neutron stars (very high density 
and low temperature), supernova explosions (high den- 
sity and moderate temperature) and relativistic heavy 
ions collisions (high temperature and not very high den- 
sity). 

It is possible, at least in principle, to obtain informa- 
tion about the behavior of matter under extreme con- 
ditions through lattice QCD calculations. Recently 
it was shown that the deconfinement transition seen at 
zero baryon density becomes a smooth crossover at very 
small density and that at low enough temperature chi- 
ral symmetry remains broken at all densities. Although 
these results need to be confirmed by other calculations, 
it is anyway interesting to explore a model where these 
features are implemented. In particular, if a constituent 
quark model is used, the deconfinement transition need 
not to be as discontinuous as in MIT-like models, where 
the transition to quark matter implies also chiral sym- 
metry restoration. 

In this letter we compute a finite temperature EOS 
based on the transition to quark matter, and study in 
particular its connection to the problem of supernova ex- 
plosion. The most recent calculations show that using 
traditional EOS a successful explosion is not achieved 
via the prompt shock mechanism and late neutrino trans- 
port still have to be more exhaustively investigated 
The prompt explosion mechanism needs a very soft EOS, 
what seems apparently incompatible with the observed 
masses of neutron stars. The last conclusion is based 
essentially on the results coming from the phenomeno- 
logical BCK EOS §,§. This difficulty can be circum- 
vented. In order to obtain a successful explosion one 
needs a softening of the EOS at densities slightly larger 
than the nuclear matter saturation density p |j>|. On 
the other hand, the maximum mass of a neutron star de- 
pends on the very large density behavior of the EOS and 
requires a not too soft EOS in that density region. 

The idea to correlate the softening of the EOS to the 



presence of a phase transition [pi and in particular to 
the formation of quark matter |7|] is rather old. This 
possibility, in connection with the problem of supernova 
explosion, has been considered recently by Gentile et al. 
H in a phenomenological way, without any attempt to 
relate the parameters governing the transition to other 
quark model calculations. Moreover the dependence of 
the EOS on the temperature and on the electron fraction 
was not investigated. Since during the first seconds of the 
its life a proto-neutron star deleptonizes, it is particularly 
important to study how the transition to quark matter 
is affected by the electron fraction. 

We present here a calculation of the transition to quark 
matter at finite temperature. In order to satisfy Gibbs 
conditions when more than one charge is conserved, one 
has to use the technique developed by Glendenning [||. 
Until now, the only calculation incorporating this tech- 
nique at finite temperature was the study of the liquid- 
gas phase transition in nuclear matter made by Miiller 
and Serot Jh}. 

We will use the Color Dielectric Model (CDM) 
to describe quark matter. The CDM is a confinement 
model which has been used with success to study prop- 
erties of single nucleons, such as structure functions ]l4|] 
and form factors Jig ], or to investigate zero temperature 
quark matter []l6f and the structure of neutron stars [[17) . 

An important feature of the CDM is that effective 
quark masses are of the order of 100 MeV even at very 
large densities, hence chiral symmetry is broken and the 
Goldstone bosons are relevant degrees of freedom. This 
is to be contrasted with models like the MIT bag, where 
quarks have masses of a few MeV. 

The Lagrangian of the model reads: 

C = i^d^ + i(o» 2 + ±(d^) 2 - U (a, n) 
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where U (a, ir) is the "mexican-hat" potential, as in Ref. 
p5| . The Lagrangian C describes a system of interacting 
u, d and s quarks, pions, sigmas and a scalar-isoscalar 
chiral singlet field x- The scalar field x characterizing the 
CDM is related to the fluctuation of the gluon condensate 
around its vacuum expectation value In the model 
this fluctuation is rather small indicating a smooth tran- 
sition between the exterior and the interior of the nucleon 
and allowing a soft deconfinement. 
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The coupling constants are given by g Ut d — g{f n ± £3) 
and g s = g(2f K - f n ), where f n = 93 MeV and fx = 113 
MeV are the pion and the kaon decay constants, respec- 
tively, and £3 = fK ± ~fK° — —0.75 MeV. These coupling 
constants depend on a single parameter g. 

Confinement is obtained via the effective quark masses 
m u ,d = 9u,dd- 1 '(x/tt) and m s = g s /x, which diverge out- 
side the nuclcon. Indeed, the classical fields x and o are 
solutions of the Euler-Lagrange equations and x S oes 
asymptotically to zero at large distances. 

In the following we will use for the model parameters 
the values: 

.9 = 0.023 GeV, M = 1.7 GeV, (2) 

giving a nucleon isoscalar radius of 0.80 fm (exp.val.=0.79 
fm) and an average delta-nucleon mass of 1.129 GeV 
(exp.val. =1.085 GeV). A similar set of parameters has 
been used to compute structure functions and form 
factors pl| and to study neutron stars |l7j . 

We describe the hadronic phase with a relativistic field 
theoretic model of the Walecka type [jl8|], including pro- 
tons and neutrons only. The parameters used to define 
the Lagrangian of the hadronic part are the ones labeled 
HS81 in the work by Knorren et al. jlj|. 

The transition to quark matter has been studied using 
the technique developed by Glendenning J9| , since in the 
transition two quantities are conserved, the baryon (B) 
and the electric (C) charge : 

PB = (l-x)PB+Xp q B (3) 
PC = (1 - X)Pc + XPc + Pe+ Pn=0. 

Here x is the fraction of matter in the quark phase. The 
superscripts h and q label the density in the hadronic 
and in the quark phase, respectively. The electron (p e ) 
and the muon (p n ) charge densities contribute to make 
the total electric charge equal to zero. Due to the pres- 
ence of strange quarks electron and muon densities are 
suppressed. 

Since matter has to be in chemical equilibrium under 
/3-decay and deconfinement, the following equations have 
to be satisfied: 

pn pp pe Pi*e : Pn Pp pn Pv^ 
%Pd + Pu = Pn , Pu — Pd = Pp — Pn 

Ps= Pd, (4) 

together with the usual condition for mechanical equilib- 
rium, i.e. the equality of the pressure in the two phases: 

P h = P q . (5) 

Other conditions depend on the specific problem under 
discussion. For instance, in order to study the structure 
of the star after deleptonization, one assumes that neutri- 
nos can escape freely, so their chemical potential is set to 
zero in previous equations. This assumption is incorrect 



in the first seconds of the life of the proto-neutron star. 
Due to neutrino opacity lepton numbers are conserved, 
neutrino chemical potentials are different from zero and 
the EOS is computed for fixed values of the lepton frac- 
tions: 

Yi e = (p e +p„ e )/p B , Y tfl = (p fi + p Vtl )/p B • (6) 

Since before the collapse Yi = 0, this quantity has been 
kept fixed and the EOS has been computed for various 
values of Yi e . 

In the computation of the EOS all previous equations 
have been solved together with the mean-field equations 
of the Walecka model and of the CDM. Finite temper- 
ature has been taken into account using the standard 
technique developed e.g. in Ref. |2(| . 
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FIG. 1. Boundaries separating hadronic matter, mixed 
phase and quark matter in the density-temperature plane. 
The labels indicate various values of Yi c , s is for symmetric 
matter. 

We come now to the results. In Fig.|] we present the 
boundaries separating hadronic matter from mixed phase 
and the latter from pure quark matter. The labels corre- 
spond to various values of Yj e . Symmetric nuclear mat- 
ter is also presented. The transition region depends on 
the electron fraction Yi e . In symmetric matter at low 
temperatures the mixed phase forms at p = 0.23 fm -3 , 
therefore no quark matter is present in heavy nuclei. De- 
creasing the value of Yi c the phase transition starts at 
lower densities. At any value of Yi e the mixed phase ex- 
tends on a rather limited range of densities and even at 
zero temperature pure quark matter phase is reached at 
densities slightly larger than 2po- At higher tempera- 
tures the transition starts at lower densities. We have 
explored only the region of moderate temperatures, the 
relevant one for supernova explosion. In order to study 
higher temperatures it would have been essential to take 
into account quantum corrections beyond mean field ap- 
proximation. 

The scenario depicted in Fig.[l] is rather different from 
the one described by the MIT model, where the transi- 
tion, at least in symmetric nuclear matter, starts at very 
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large densities and the mixed phase extends on a broad 
density range. Using the MIT model, which assumes the 
interior of the nucleon to be in a perturbative regime, it 
is almost unavoidable to have a first order transition. In 
the CDM it is conceivable, at least in principle, to ob- 
tain a smoother transition by taking into account quark 
correlations beyond mean field approximation. 
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TABLE I. Total binding energy, gravitational and internal 
energies in Walecka model and in CDM (see text). All energies 



are in foe (10 erg 




FIG. 2. Pressure (upper box) and adiabatic index (lower 
box) as function of the density in CDM (solid) and in BCK 
(dotted). The pressure in Walecka model is also shown 
(dashed). 

To investigate our EOS in connection with the prob- 
lem of supernova explosion, we compare with BCK EOS 
H (parameters as in model 38). The latter is a totally 
phenomenological EOS which is soft enough to allow for 
supernova explosion, but gives a maximum mass smaller 
than the mass of PSR 1913+16 (1.44 M ) |§. In our 
model the maximum value of the gravitational mass for 
a non-rotating cold star is 1.59 M Q Jl7| . 

In Fig.|| we present results for Yi e = 0.4 and entropy 
per baryon number S/R = 1. Due to the presence of 
strange quarks, the electron fraction Y e is rather small 
in the quark phase, Y e ~ 0.3, and the muon fraction 
is always very small. In the upper box we compare the 
pressure in the Walecka model, in our model and in BCK 
EOS. Due to the phase transition, our EOS is rather soft 
from p = 0.17 fm -3 to p = 0.34 fm -3 . On the other 
hand, after p = 0.34 fm -3 it is considerably stiffer than 
BCK, allowing higher masses for the proto-neutron star. 
These conclusions are strengthened by the computation 
of the adiabatic index, shown in the lower box of FigJ|. 
Clearly in the mixed phase matter offers little resistance 
to collapse, but when pure quark matter phase is reached 
the collapse is halted. In the mixed phase region our 
adiabatic index is even smaller than in BCK. 

Another way of investigating the possibility of a 
prompt explosion using our EOS is to estimate the 
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extra energy at disposal for the shock wave in our model 
respect to pure Walecka model. To this purpose we follow 
the analysis of Gentile et al. @]. Their way of estimat- 
ing the extra energy is the following: they compute the 
binding energy of the central part of the proto-neutron 
star, immediately before bounce, integrating on a mass 
of 0.5Mq. The binding energy is the sum of gravita- 
tional and internal potential energy: BE = V g + Ei = 
gravitational mass - baryonic mass. Integrating only on 
the inner region, a positive value for the BE is obtained, 
indicating that this region tries to expand, pushing the 
exterior envelope. The larger the (anti-)binding energy, 
the stronger the push. We can therefore compare these 
numbers in Walecka and in our model, using the EOS 
with S/R — 1 and Yj e = 0.4. The result is shown in the 
Table. An extra energy of about 8 foe is at disposal for 
the explosion. 
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3. Composition of a star having S/R = 1 and 



Yi c = 0.4. A different length's scale has been used where 
the mixed phase is formed. 
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FIG. 4. Same as FigJ^, for a star having S = and /3-stable. 

We come now to the second problem we like to study, 
namely the structure and evolution of the proto-neutron 
star. In FigjH we show the composition of a star with 



Yi c = 0.4 (Y e ~ 0.3) and entropy per particle S/R = 1. 
The central temperature is of the order of 10 MeV. These 
conditions should be realized at the bounce. Immediately 
after the bounce the entropy increases. After a time of 
the order of 10 seconds the proto-neutron star cools down 
and deleptonizes. After that the composition of the star 
does not change any more. We show this later stage in 
Fig.||, where we assume entropy S — and /3-stability. 
In both figures the baryonic mass is the same, 1.54M Q , 
which corresponds to a gravitational mass of 1.4M0 for 
the final star of Fig.[|. 
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FIG. 5. Mass fractions in hadronic, mixed and quark phase 
(upper box) and binding energy (lower box) as function of Yi e 
for entropy S/R = 1 (solid) and S/R = 3 (dashed), respec- 
tively. 

There have been many speculations about the possi- 
bility of a late neutrino emission, based on the SN1987 
data pljj . In particular the idea of a late transition to 
strange matter with a new emission of energetic neutri- 
nos has been invoked. We investigate this possibility in 
Fig.||, where we show the composition of the star and its 
binding energy as a function of the lepton fraction Yi c 
and for two values of the entropy per particle. To study 
this problem one should really solve the dynamics' equa- 
tions giving lepton fraction and entropy as a function of 
time. What can be learnt from Fig.[| is the absence of 
a sudden jump in the composition of the star and in its 
binding energy as a function of Yi e . The conclusion we 
can draw is that an emission of energetic neutrinos dur- 
ing the first seconds is indeed possible, since the binding 
energy increases steadily as Yi c is diminishing. Anyway 
a peak in the neutrino luminosity can be due only to the 
dynamics of the explosion, not to a discontinuity in the 
EOS of the proto-neutron star. 

Recently there has been a speculation about the pos- 
sibility of having a proto-neutron star collapsing later 
to a black hole, due to the softening of the EOS after 
deleptonization In our calculation we do not have 

this effect, because the maximum baryonic mass for the 
proto-neutron star is essentially independent of the value 
of Y e and S and remains always near 1.82M . 
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In this letter we have shown the possibility to have 
a phase transition at densities slightly larger than po in 
pre-supernova matter. This transition softens the EOS 
which in this range of densities is comparable with the 
phenomenological BCK EOS, but gives acceptable values 
for the maximum gravitational mass of the final neutron 
star. Our EOS indicates also the possibility of a energetic 
neutrino emission in the first seconds of the life of the 
proto-neutron star. All these conclusions will soon be 
tested in a dynamical simulation of the explosion. 

It is a pleasure to thank V.Barone, Z.Berezhiani, 
L.Caneschi, H.-T. Janka and W.Keil for many stimulating 
discussions. 
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